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Abstract. Rates of convergence of solutions of various two-dimensional a— regularization models, sub- 
ject to periodic boundary conditions, toward solutions of the exact Navier-Stokes equations are given in 
the L°°-L^ time-space norm, in terms of the regularization parameter a, when a approaches zero. Fur- 
thermore, as a paradigm, error estimates for the Galerkin approximation of the exact two-dimensional 
Leray-a model are also presented in the L°°-L^ time-space norm. Simply by the triangle inequality, one 
can reach the error estimates of the solutions of Galerkin approximation of the a-regularization models 
toward the exact solutions of the Navier-Stokes equations in the two-dimensional periodic boundary 
conditions case. 
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1. Introduction 

The closure problem of averaged quantities in turbulent flows has been, for many years, an out- 
standing open challenge for turbulence models. In the recent decade, various a-regularization models 
(Navier-Stokes-a, Leray-a, Modifled Leray-a, Clark-a and simplifled Bardina model) were introduced as 
efficient subgrid scale models of the Navier-Stokes equations (NSE) [2, 6, 10, 11, 12, 15, 18, 19, 25, 26, 27] 
(see also [37] for an analytical study of a mathematical generalization of the Navier-Stokes— a model). 
In particular, it was shown in some of these papers that these a-models fit remarkably well with em- 
pirical experimental data for a large range of huge Reynolds numbers. Moreover, these models were 
implemented numerically by various groups [3, 4, 12, 13, 20, 21, 22, 24, 25, 28, 31, 32, 33, 36]. Indeed, 
the authors of [30] have pushed further this numerical analysis point of view, concerning the a-models, 
in their study of the MHD-a model (sec also [34, 35]). In fact, there have been extensive analytical 
studies on the global regularity of solutions and finite-dimensionality of global attractor of these models, 
however, there is much less work on the convergence, especially, the rate of convergence of solutions of 
various a-models toward the solutions of the exact NSE, when the regularization parameter a approaches 
zero. The authors of [19] study the convergence of the three-dimensional (3D) Navier-Stokes-a (NS-a) 
model to the 3D NSE. To be more specific, they show that there exists a subsequence of solutions m"^ 
of the 3D NS-a that converges to one of the Leray-Hopf weak solutions of the 3D NSE, with periodic 
boundary conditions. Similar results are reported in [30] concerning the MHD-a model. Later the au- 
thors of [14, 45] show that the trajectory attractor of the Leray-a and Navier-Stokes-a, respectively, 
converges to the trajectory attractor of the 3D Navier-Stokes system, as a approaches zero. Since the 
uniqueness theorem for global weak solutions (or global existence of strong solutions) of the 3D NSE is 
not yet proved, the studies mentioned above either consider convergence to a weak solution or consider 



Date: October 14, 2009. 

To appear in: Numerical Functional Analysis and Optimization. 

1 



2 



Y.CAO AND E.S. TITI 



the convergence to the trajectory attractor of weak solutions. Recently the authors of [9] study the con- 
vergence rate of the Navier-Stokes-a model and obtain a mixed L^-L^ time-space norm for small initial 
data in Besov-type function spaces in which global existence and uniqueness of solutions for 3D NSE can 
be established with "small" initial data and external forcing. Similar results can also be derived by ap- 
plying the same techniques for the other a-models under the assumption of existence of strong solutions 
of the 3D NSE. e.g.. when the initial value and external forcing arc small enough in the appropriate norms. 

It is worth mentioning that inspired by the a-regularization models of turbulence similar regulariza- 
tion schemes were introduced and implemented in [7, 8], in the context of the nonlinear Schrodinger 
equation, and in [4] for the Leray-a regularization of the inviscid Burgers equation. 

We mainly investigate in this paper the rates of convergence of four a-models (NS-a model, Leray-a 
model. Modified Leray-a model and simplified Bardina model) in the two-dimensional (2D) case, subject 
to periodic boundary conditions on the periodic box [0,L]^. Since unique strong solution is proved to 
exist globally in time with any smooth enough initial data in the 2D case, we will show upper bounds, 
in terms of a, for the difference between solutions of the 2D a-modcls, u", and solutions of the 2D NSE 
system, u, in the L^-norm for any time interval [0,T]. Specifically, we show that all the four a-modcls 
we include in this study have the same order of convergence and error estimates, i.e., the L^-norms of 
the differences, \\u — u"||2, are of the order O (( j-)(log(^))^/^) , as ^ tends to zero. These results are 
presented in detail in Section 4. It is worth mentioning that the Brezis-Gallouet inequality plays an 
essential role in our error estimates in the 2D case; which in turn results in the logarithmic factor. This 
logarithmic factor, however, is absent in the 3D case which is treated in [9]. In a forthcoming paper, 
we will consider the rate of convergence of 3D Leray-a and NS-a models toward the adequate strong 
solution of the 3D exact NSE system, provided the latter exists, and compare the results to that of [9]. 

In Section 5, we consider the error estimates of the Galerkin approximation solutions in the 2D case, i.e., 
we estimate the difference between the solutions of the 2D a- model, , and solutions of its corresponding 
finite-dimensional Galerkin approximation system, w"^, where m > 1 is the order of the truncation mode 
(dimension) of the Galerkin system. We will study, as an example, the Leray-a model and present the 
detailed proof of the error estimates. One can easily apply similar arguments to the other a-regularization 
models (NS-a model. Modified Leray-a model and simplified Bardina model) and obtain similar error 
estimates for these models. For the Leray-a model, the i^-norm of the difference, ||it" — Wm||2, is of 

the order O ^^|-p-(log(Am+ii^))^/^j , under the assumption that a is small such that a^ < ^ , 

where A^-i-i is the (m + 1) — th eigenvalue of the Stokes operator in the 2D case. Applying the triangle 
inequality, we get error estimates concerning solutions of finite-dimensional Galerkin system of the 2D 
Leray-a model, w^, as an approximation of the exact solution of the 2D NSE system, u. Specifically, 
we show that the error in the i^-norm of the difference between the solution u of the 2D NSE and m^, 

the solution of the Leray-a Galerkin system, ||m — u^||2, is of the order O (^j^^-^{\og{Xm+iL^)^^'^)j , 

provided a is small enough satisfying a < j-^^^- 

Before we present the main results of the error estimates, we will first introduce in Section 2 some 
notations and preliminaries that will be used throughout this paper. In Section 3 we will present all four 
a-models in functional setting and establish a priori estimates for the solutions, by investigating the 
finite-dimensional Galerkin systems then passing to the limit to — ^ oo by Aubin compactness theorem to 
obtain upper bounds for the exact solutions of the relevant models in certain norms. The main results 
for the rate of convergence, in terms of the regularization parameter a, of the a-modcls in the 2D cases 
will be presented in Section 4. After establishing the rate of convergence of solutions of a-regularization 
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models toward the solutions of the exact Navier-Stokes system, we will further show in section 5 the er- 
ror estimates of difference between solutions of the 2D Leray-a model and solutions of its corresponding 
finite-dimensional Galcrkin approximation system. 

2. Functional Setting and Preliminaries 

In this section, we introduce some preliminary background material following the usual notation used 
in the context of the mathematical theory of Navier-Stokes equations (NSE) (see, e.g., [16], [40], [41], 
[42]). 

(i) We denote by L'^ and i?'" the usual Lebesgue and Sobolev spaces, respectively. And we denote 
by I • I and (■, ■) the L^-norm and L^-inner product, respectively. 

(ii) Let J- be the set of trigonometric polynomials of two variables with basic periodic domain 

— [0, LY and spatial average zero, i.e., for every (j) G (j){x) dx = 0. We then set 

V= [Tf : V-0 = O}. 

We denote by H and V the closures of V in the [L^Y ^^'^ (^^)^ topologies, respectively. For 
u E H and v E V , we denote by 

\u\^ = I \u{x)\^ dx and = / \Vv{x)\^ dx 

Jn Jn 

the norms in H and V, respectively. We also note that by Rellich Lemma (see, e.g., [1]), we 
have that V is compactly embedded in H. 

(iii) For any s > 0, we denote by if^ the closure of J- with respect to the H^{fl) topology. Hence, 
for any u G , we can write the Fourier expansion 

u{x) = E Tifce'^^"^, 

where the Fourier coefficients iik satisfy the reality condition ii^k = w^, Vfc G Zq := Z^\{0,0}. 
We define the norm on this space as 

iK. = (^f<-'-^\s (1 + 1^1^)1^,1^ 

(iv) We denote hy Pa : ^ H the Helmholtz-Leray orthogonal projection operator, and hy A = 
—P^A the Stokes operator, subject to periodic boundary conditions, with domain D{A) = 
{H'^{n)Y n V . We note that in the space-periodic case 

= -P^Au = -Au, for all u G i:'(A). 

The operator is a self-adjoint positive definite compact operator from H into H (see, e.g., 

[16], [41]). We denote by < (^f = M < M < the eigenvalues of A in the 2D case, 

repeated according to their multiplicities. It is well known that in two dimensions the eigenvalues 
of the operator A satisfy the Weyl's type formula (see e.g., [16], [42]), namely, there exists a 
dimensionless constant cq > such that 

-<^<coj, forj = l,2,... (1) 

Co Ai 

We also observe that in the periodic case, D(A"/2) = (iJ"(Sl))2 n H, for n > 0. In particular, 
one can show that V = D{A^^^) (see, e.g., [16], [41]). 

(v) For every w E V, we have the Poincare inequality 

Ai|^|^<||«;f. (2) 
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Moreover, one can easily show that there is a dimcnsionless constant c > 0, such that 

c\Aw\ < \\w\\h2 < c~'^\Aw\ for every w G D(A), (3) 

and, by virtue of Poincarc incquahty, 

c|A^/^w| < \\w\\hi < c^^\A^^^w\ for every w &V. (4) 

Hereafter, c wiU always denote a generic dimcnsionless constant. 

Notice that, thanks to (4). the norm of V is equivalent to the usual norm. 

(vi) Let {wj}J^i be an orthonormal basis of H consisting of cigenfunctions of the operator A. Denote 
by Hm = spanjwi, W2, ■■■,Wm,}, for m > 1 and let Pm be the i^-orthogonal projection from H 
onto H„i, then it is easy to see that 

\{I-P„,)cf,\-' <-^\\cf,\\^, foraU^GF, (5) 

Am+l 

\\{I - P^)cj)\\^ < -^\A(j>\\ iov all (f> e D{A). (6) 

Moreover, one can also easily show that 

\APmcj)\^ < A,„||F,„0f , for all G D{A). (7) 

(vii) We recall the following 2D interpolation and Gagliardo-Nirenberg-Ladyzhenskaya inequality 
(see, e.g., [1], [16], [41]) : 

ll0lU^<c||</>|li{'||0||^/f. (8) 

(viii) For every Wi,W2 G V, we define the bilinear operators 

B{WUW2)^PA{W1-S/)W2), (9) 
B{wi,W2)^~PAwiX{y XW2)). (10) 

In the following lemma, we will list certain relevant inequalities and properties of B (see, e.g., [16], [41], 
[42]) and of B (see [19]). 

Lemma 1. The bilinear operator B defined in (9) satisfies the following: 

(i) B can he extended as a continuous bilinear map B : V x V V' , where V' is the dual space of 
V. 

In particular, the bilinear operator B satisfies the following inequalities: 

\{B{u,v),w)v' \ < c\u\^ \\u\\^ \\v\\\w\^ \\w\\^ , for all u,v,w G V, (11) 

\{B{u,v),w)\<c\\u\\oo\\v\\\w\, for allu E D{A), V eV, w e H, (12) 

\{B{u,v),w)\<c\u\\\yv\\^\wl for allu e H, V e D{A^/'^), w e H, (13) 

\{B{u,v),w)(^DiA)y\<c\u\\\v\\\\w\\^, forallueH,veV,weD{A). (14) 
Moreover, for every wi,W2,W3 G V, we have 

{B{wi,W2),W3)v' = -{B{wi,W2),W2)v', (15) 

and in particular, 

{B{wuW2),W2)v' ^0- (16) 

(ii) In the 2D periodic boundary condition case, we have 

(i?(^,^),A^) = 0, (17) 

for every ip G D{A). 
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(iii) B can also be extended as a continuous bilinear map B : V x V ^ V . Furthermore, for every 
Wi,W2,W3 € V, the bilinear operator B satisfies the following inequality: 

\{B{w^,W2),W3)v'\ < c\w,nw,\\'/^W2\\\w,nw3\\'^\ 
Moreover, for every wi,W2,W3 € V, we have 

{B{wi,W2),W3)v' = {B{wi,W2),W3)v' - {B{w3,W2),Wi)v', (18) 

and consequently, 

{B{wuW2),wi)v' ^0. (19) 
Lemma 2. For every u G D{A) and w ^ V , we have 

\{B{w,u),Au)\ < c\\w\\\\u\\\Au\. 
Proof. Let w eV and u e D{A). Then 

2 /■ a 

'j,™=iJo OX J dxj^ 



{B{w,u),Au) = - T. I w-j{^^ui){-^-^ui) dx 



2 

l,j,m=l \_Jq^ dXm '''"'dXj '^''dXm Jn ''^dXjdXjn dXn 

By relation (16) the second term on the right-hand side above is zero. Therefore, by the Cauchy-Schwarz 
incquahty 

\{B{w,u),Au)\<c\\w\\\\Vu\\l,, 

and by (8), (3) and (4) we have 

\{B{w,u),Au)\ < cllwllllulllAwl. 
From the above and the density of V in we conclude our lemma. □ 

Next, we state a two-dimensional periodic boundary condition version of the well-known Brezis- 
Gallouet inequality [5]. For the sake of completeness, we will present the proof of this version in the 
Appendix. 

Proposition 3. There exists a scale invariant constant c > such that for every ip £ Di^A), 

1/2 



Applying the Brezis-Gallouet inequality above to (12) and (14), we have the following corollary (see also 
[43] for similar logarithmic inequalities concerning the bilinear term). 

Corollary 4. In the two-dimensional case, the bilinear operator B satisfies the following inequalities: 
\{B{u,v),w)\<c\\u\\\\v\\\w\(l + log-^^^] for allueD{A),veV,weH, (20) 



,|2 \ 1/2 

where Ai = (x")^ /irsi eigenvalue of the Stokes operator A 



I {B{u,v),w) I < c|u|||w|l||u;|| ( 1 + log ) , for all u e H, v e V, w e D{A). (21) 



3. A Priori Estimates of the 2D Navier-Stokes Equations and the a-REGULARizATiON 

Models 



In this section, we will establish a priori estimates for solutions of the 2D NSE and the 2D a- 
regularization models. These results will be useful for the error estimates in the next two sections. 
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3.1. A priori estimates for the 2D Navier- Stokes equations. We recall that the two dimensional 
(2D) NSE is equivalent (see, e.g., [16], [41]) to the functional evolution equation in the Hilbert space H 

nil 

— + vAu + B{u,u) = /, (22) 
u(0) = uq. 

The corresponding Galerkin system of the 2D NSE is given below as a system of ordinary differential 
equations in the space Hm, defined in Section 2: 

— + vAu,n + PmB{u ) = fm, (23) 

Um{0) = Uom- 

where Wom = Pm "O, fm = Pm f ■ 

The proof of existence and uniqueness of solutions of the above 2D NSE system (22), subject to periodic 
boundary conditions, can be established by applying Galerkin approximation procedure [16, 41, 42]. The 
idea is to establish a priori estimates for the solutions of the finite-dimensional system (23) and then by 
applying Aubin compactness theorem (see also [16, 41, 42]) one can extract a subsequence that converges 
to the unique solution of the NSE system (22). The details of the proof are textbook material and will 
not be presented here. In the following proposition we will establish a priori estimates for solutions of 
the Galerkin system (23) that we will use later. 

Proposition 5. Let uq G H,T > and f G L^{{0,T);H). Let Um he the solution of the Galerkin 
approximation of the 2D NSE, system (23), over the interval [0,T] with initial data uom for a given 
m > 1, then 

sup (Mt)f + uf \\u,n{s)f ds] <K^{iy,uoJ,L,T), (24) 

0<t<T \ Jq / 

where |uoP + ^ |/(s)P ds and Ai = {^f. 

Proof. We take the inner product of the first equation in (23) with to obtain 

]^^^\Um? + l^\\u,nf = iPmf,U„,) < \ f\\u^\. 



By Poincare inequality (2) we have 



Hence, 



1^1 II 112^ l/lll"m|| ^ \f\' ^^1, 



^\u„,\^ + iy\\u„,f < 

dt v\i 



and by integration we have 

I Urn (OP 

which concludes the proof. □ 



vl ||M,„(s)f ds< |uo|' + ^ / |/(s)|'ds foralHe [0,r], 



Remark 1. One needs to further establish estimates on and then, thanks to Aubin compactness 
theorem [16, 41], one can extract a subsequence Mm', that converges to the unique solution u of system 
(22), as m' oo. Following the standard procedure as in [16, 41], one can show that the solution of the 
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NSE (22) satisfies the same a priori estimates, namely, let u be the solution of the system (22) over the 
interval [0, T] with initial data uo G H, then 

sup (\uit)\^ + iyf \\uis)\\^ds) <K^iiy,uoJ,L,T). (25) 

0<t<T \ Jo / 

3.2. A priori estimates for the 2D Leray-a model. In this subsection,we will establish a priori 
estimates for solutions of the 2D Lcray-a rcgularization model. The Leray-a model was introduced and 
analyzed in [15] and implemented computationally [15, 20, 21, 22, 31] in the context of subgrid scale 
models of 3D turbulence (see also [4, 24] and [31] for 2D computations with the Leray-a model). The 
Lcray-a model was inspired by the Navicr-Stokcs-a model (also known as viscous Camassa-Holm or 
Lagrangian-averaged-Navier-Stokes-a model) [10, 11, 12, 18, 19], and it happened to fit as a member 
of the general family of regularizations introduced in the seminal work of Leray [29] in the context of 
establishing the existence of solutions for the 2D and 3D NSE. 

The Leray-a rcgularization model of the NSE is given by the following functional evolution system 
in the space H: 

— +vAv^+B{u",v^) = /, (26) 
u"(0) = ua. 

Observe that when the rcgularization parameter a = one recovers the exact NSE system (22). 

The corresponding Galerkin system of the 2D Lcray-a model is given below as a system of ordinary 
differential equations in the space iJ,™: 



dt 



i.v4«;^ + P™i?«„0==/™, (27) 

"m(0) = "Or, 



m 1 



where Uom = Pm "O, /m = Pm f- 



The proof of existence and uniqueness of solution of the above Leray-a system (26) and the other three 
a-models we cover in this section, subject to periodic boundary conditions, can be established by apply- 
ing Galerkin approximation procedure. One can follow similar steps as those for the NSE [16, 41, 42] to 
show the proof. The proof is not the heart of this paper and we will omit the details and only establish 
a priori estimates for solutions of the Galerkin system (27). To be concise and focus on the essential 
matter of this paper, in this and the following subsections of a priori estimates for the 2D a-models, we 
will simply skip the details of the proof and will not restate these comments. Interested readers can refer 
to the relevant literature in [16, 41, 42] and references therein to fill in the gap. 

Proposition 6. Let uo £ D{A), T > and f e L^((0, T);H). Let he the solution of the system (27) 
over the interval [0,T] with initial data ito,„ for a given m>l, then 



sup 

0<t<T 



|C(Of + «'l^<.(i)P + ^ / V<„(s)P + aM"^<{s)\^)ds 

Jo 

where ||^io||' + a'|A^.o|' + ^ Jo 1/(^)1' ds. 
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Proof. We take the inner product of the first equation in system (27) with Au"^ to obtain 



By (16) and (17) one has 

{B{K„vZ),Au^„) = (i?«„<J,A<J +a2(i?«,,A<J, A<J = 0. 
Consequently, from the above we have 



< 



\fl 

2iy 



-\Au 



a |2 



Hence 



and by integration over the interval (0, t), for t G [0, T] we have 



||<„(t)||^ + a1^<„(t)|^ + W (|^<(.)|^ + a^|^^/X.(^)r)rf^ 

Jo 

|2 j„ , 11 , 112 , „,2| 



< 



1 



\f{s)\'ds+\\uo„ 



a \Auor, 







1 

< - 



\f{s)\^ds+\\uof + a^\Auo\, 



which concludes our proof. 



□ 



Remark 2. Similar to the 2D NSE, one needs to further establish estimates on and then, thanks to 
Aubin compactness theorem [16, 41], one can extract a subsequence u'^, that converges to the unique 
solution of system (26), as m' —> oo. As a result, one can also prove that the solution of the Leray-a 
system (26) satisfies the same a priori estimates as in Proposition 6, namely, let be the solution of 
system (26) over the interval [0,T] with initial data uq £ D{A), then 



sup 

0<t<T 



\\u"{t)\\^ + a''\Au"it)\^ + iy / i\Au''is)\^ + a^\A^^^u"{s)\^)ds 



(28) 



We emphasize here again that the details of the proof of the above results are omitted. In the following 
subsections of a priori estimates of the NS-a, Modified Leray-a and simplified Bardina models, we will 
also skip the details and will not restate the remark. 

3.3. A priori estimates for the 2D NS-a model. In this subsection, we will establish a priori 
estimates for the 2D NS-a regularization model of the 2D NSE. The NS-a model (also known as the 
viscous Camassa-Holm or Lagragian-averaged-Navier-Stokes- a model) was introduced and analyzed in 
[10, 11, 12, 18, 19, 25], which was also the first of the family of the a-models. In addition to the remark- 
able match of explicit solutions of the 3D NS-a model to experimental data, in the channels and pipes, 
for a wide range of huge Reynold numbers [10, 11, 12], the validity of NS-a model as a subgrid scale 
turbulence model was tested numerically in [12, 13, 22, 24, 32, 33, 36]. 



The NS-a regularization model of the NSE is given by the following functional evolution system in 
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the space H: 



+ iyAv°' + B{u°',v°') 



(29) 



ii"(0) = uo, 

where B is given by (10) and Lemma 1. 

Observe that when the regularization parameter a — one recovers the exact NSE system (22). 

The corresponding Galerkin approximation system of the 2D NS-a model is given below as a system 
of ordinary differential equations in the space Hm- 



<n(0) 



fm, 



(30) 



where WOm = Pm Uq, /,„ ^ Pm f . 



The proof of existence and uniqueness of solutions of the NS-a system (29), subject to periodic boundary 
conditions, is similar to that of the 2>D case as it is presented in [19]. In the next proposition we will 
omit the details and only show the a priori estimates for solutions of system (30). 

Proposition 7. Let uq € D{A),T > and f G L^{{Q,T); H). Let uj^ be the solution of the Galerkin 
system of 2D NS-a (SO) over the interval [0,T] with initial data uq„i for a given m > 1, then 



sup 

0<t<T 



<.(op+«'ii«™wip+^ / {\K,{s)r+aM<{sr)ds 



(31) 



where Kqj^ :— |uoP + q;^||uo|P + 7^ /q |/(s)P ds, which is independent of m, and Ai = (^)^ 
Moreover, 



sup 

0<t<T 



\mt)r + a'\AKM' + '^ / {\AKA^)\'+a^\A'/^u'^{s)\')ds 



< K 



02 1 



(32) 



where ll'"o|P + Q:^|^uoP + j; Jq \fis)\'^ ds + -^KqqKqi, and Kqq is « constant depending on 

Aitol,!^, r, / and Kqi, which is give explicitly in (36). 

Proof. By taking the inner product of the first equation in (30) with and using (19) we obtain 



1 d 
2dt 



{\u 



a |2 



ai!<r)+Kii<.r+«i^<„i^) 



iPmf,u"J<\f\\u"J< 



i/iii«™ii 



1/2 



< 



lit 

2v\i 



,a 112 



+ 



where in the last two steps we used Poincare inequality (2) and Young's inequality. 
As a result we have 

Integrating the above in time we obtain (31). 
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Next, we prove (32). We take the inner product of the equation in (30) with Au'^ and follow simi- 
lar steps as above to obtain 

1 d 



2dt 



iB^f^Au'^J<l±^ + -\Au'^\^. 



Observe that 

where we use in the above (18), (15), (16), (17) and the relation v"^ ~ + a^Au'^. 
From (34) and Lemma 2 we have 



|(B(u^,«,:J.),^z.^)| <c«2|A3/2^.^||i<Jl|^<„| < -a'lA'/'u'^J' + —\\u-J'\A<\'- 
From the above and (33) we obtain 

|(l|u^f + a2|Au^l2) + z.(|A<„p + a2|^3/2^a|2) 



< 



I/I 



ca 



+ — iKJii^cr 

- ^ + ^(IKf + «'l^"™l')'- 



Therefore, by Gronwall's inequality we have 

\\u^it)\\'+a'\Au^it){' < e^/o(ll"™(^)lP+"^l^"™(^)l^)'^«(||^o||2 + „2|^^^|2) 



+ i / e^/i("""('')ll'+"'l^"™('')l')''^|/(a)|2da, 



and by (31) we obtain 



||<,(t)r + a^Kn"(t)l' <e^^^-(||iior + alAi.on + -e^"- / |/(a)|^ da ifo'o- (36) 

^ Jo 

Integrate (35) and using (31) and (36), we obtan 

\\KMr + o^'\AK^it)\' + f\\AK,is)\^+a'\A'/^uf^{s)\')ds 



< \\uor + a'\Auo\' + - I \f\'dt + -KioKL 



c ~, 



V 



which concludes the proof. 



(33) 

(34) 



(35) 



□ 



Remark 3. By following similar steps as those of NSE [16, 41, 42], one can show that the exact solutions 
of the 2D NS-a system have the same a priori estimates as those of the solutions of the Galerkin system, 
namely, let be the solution of the system (29) over the interval [0, T] with initial data G D{A), 
then 



sup 

0<t<T 



sup 

0<t<T 







iu"(or + «iiw"(t)ir + / (h"(s)r + a^iAM"(s)i^)ds 

'0 



|lM"(t)f + a2|AM"(t)p + jy / (|Au"(s)p + a2|A3/2u"(s)|2)rfs 



02- 



(37) 
(38) 
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3.4. A priori estimates for the 2D Modified Leray-a model. In this subsection, wc will establish 
a priori estimates for the 2D Modified Leray-a (ML-a) regularization model of the 2D NSE. Inspired by 
the NS-a and Leray-a models, the ML-a was introduced and analyzed in [26] and is tested numerically 
in [22] in the context of 3D sub-grid a-modcls of turbulence. 



The ML-a regularization model of the 2D NSE is given by the following functional evolution system 
in the space H: 

+ iyAv°' +B{v'',u°') = /, (39) 



dt 



Observe that when the regularization parameter a ~ one again recovers the exact 2D NSE system (22). 



The corresponding Galerkin approximation system of the 2D ML-a model is given below as a system of 
ordinary differential equations in the space Hm- 

dv" 

- H-i.^<„ + P,„i3«„«:^) = /„, (40) 



dt 



where uom Pm uq, =^ P,^ f. 



C(0) 



UOm, 



Following similar steps as in the 3D case in [26] one can establish the global existence and unique- 
ness for the system (39) in the case of 2D periodic boundary conditions. In the next proposition we will 
present a priori estimates for the solutions of the Galerkin system (40). 

Proposition 8. Let uq G D{A),T > and f G L^{{Q,T); H). Let uf^ be the solution of the Galerkin 
system of 2D ML-a (4-0) over the interval [0,T] with initial data uom for a given m>l, then 



sup 

0<t<T 



2 I 2| 

' a 



<wr+^ / {\ms)r+a'\Au^{s)nds 





<^"oi- 



Moreover, 
sup 

0<t<T 



\u^it)r + aM<it)\' + '^ / i\Au'^is)\'+a'\A'/%^{sr)ds 

JO 



< K 



02 1 



where Kqi and 7^02 o-i^e the same as in Proposition 7. 

Proof. Let us take the inner product of the equation in (40) with to obtain 

1 d 11?;" II 

2^(KI' + «'ll«™f) + KII<JP + "'l^"^P) = (^™/,<J < I./IKI < l/l^ 



< 



,a ||2 



thus we get 







(41) 
(42) 



(43) 



Integrating the above over the interval (0, t), we obtain 

KJt)\'+a'\\u:Mr + ^ [\\KMr + a'\Au:js)\')ds < f \M1 ds + \uo\' + a^uof 







lyXi 
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(44) 



which imphes (41). 

Next, we take the inner product of the equation in (40) with Au"^ to obtain 

= -{B{vZ, u^), AufJ + {P„J, A<J 
< |(i?«,z,J^),^«^)| + ^ + ^|A<p. 

Next we estimate 

\B{v'^,u:,),Au'^;,)\ = |(i3«„,0,^zi«) + a2(B(A<,0,^<J|. 
By (17) and Proposition 7 we have 

|(S(i;;^,<J,A<J| = a2|(B(A<„,<J,^zi^)| <c«2|A3/2^,^J||y,"JP^,^J 

< ^a^|A3/2^.^p + ^||u^f(a^|A<J^) 

< 'iaM'^'uzr + ^i\\<r+c^'\A<fr- 

From above we have 

|(||<JP + a'\Au:S) + KIACP + a'\A'/'u'^\') < ^ + ^(|KCf + "'I^CH', 

which is exactly the same inequality as (35) (as in the proof of Proposition 7). Following the same steps 
we obtain (42). □ 

Remark 4. Similarly, the solutions of the exact ML-a system have the same a priori estimates as those 
of the solutions of the Galerkin system, namely, let u" be the solution of the system (39) over the interval 
[0,T] with initial data uq G D{A), then 



(45) 



sup 

0<t<T 



\u''{t)\''+a^\\u"{tW+iy / {\\u''{s)W' + a^\Au"{s)\^)ds 



Moreover, 



sup 

0<t<T 



u"(Of +a2|Au"(t)|2+zy / (|Aw"(s)|2+a2|A3/2u"(s)|2)d. 



< K 



02 1 



(46) 
(47) 



where Kqi and Kq2 are the same as in Proposition 7. 



3.5. A priori estimates for the 2D simplified Bardina model. In this subsection, we will establish 
a priori estimates for the 2D simplified Bardina model for the 2D NSE. The Bardina closure model of 
turbulence was introduced first by Bardina et al [2] and later simplified and studied further in [27] and 
in [6]. In particular, global well-posedness for the 37? simplified Bardina model was established in [27]. 
In [6] these results were slightly improved and the long-time behavior was investigated. Further, it was 
established in [6] the global regularity of the 3D inviscid version of the simplified Bardina, a property 
which is still out of reach for the other a-models of turbulence. 



The simplified Bardina regularization model of 2D NSE is given by the following functional evolution 
system in the space H: 

di)" 

+ i^Av°' +B{u°',u°') ^ /, (48) 



dt 



^Au" = 



w"(0) 
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Observe again that when the regularization parameter a = one recovers the exact 2D NSE system (22). 

The corresponding Galcrkin approximation system of the 2D simphficd Bardina model is given below as 
a system of ordinary differential equations in space Hm '■ 

dvl 



dt 



m 1 



a. 



(49) 



where uq,,, = P,„ uq, f,n ^ Pm f- 



In the following proposition we will establish a prior estimates for the solutions of the Galcrkin sys- 
tem (49). 

Proposition 9. Let uq G D{A), T > and f e L^{{0, T);H). Let u'^ be the solution of the system (49) 
over the interval [0,T] with initial data ito,„ for a given m>l, then 



sup 

0<t<T 



where Kq is given in Proposition 6. 

Proof. We take the inner product of the first equation in system (49) with Auf^ and using (17) to obtain 



Hence 



|(||<JP + a2|^«JJ.n+K|A^.™P+"^l^'/V„P)<^. 
dt V 



Integrating the above from to t, we have 



|K(i)f + o?\Au-M^ + V [\\Au^{s)\' + a2|^3/2<(s)n ds < - T |/(.)p ds + WuoW + a^uol 



which concludes our proof. 



□ 



Remark 5. Similar to the cases of the other a- regularization models, one can show that the solutions of 
the simplified Bardina model have the same a priori estimates as those of the Galerking system, namely, 
let u" be the solution of the system (48) over the interval [0, T] with initial data uq £ D{A), then 

nt 



sup 

0<t<T 



Wu'^it)]]^ + a^\Au"{t)f + iy / {\Au''is)f + a^lA'^/^u^is)]^) ds 
where Kq is given in Proposition 6. 



< K 



0' 



(50) 



4. Rates of Convergence of a-REGULARizATiON Models to the Navier-Stokes Equations 

We aim here to show the convergence rates of solutions of the various a- models toward the correspond- 
ing solution of the exact NSE equation when the regularization parameter a approaches zero. We will 
focus on the L^-norm of the difference between u", solution of the underlying a-model, and m, solution 
of the exact NSE. From the results concerning all four a-models included in this study, we observe that 
all the errors |m — u"\ are of the same order of 0( j-(log ^)^/^), as goes to zero. Though the four 
a-models we investigate in this paper share the same order of error estimates with the only difference 
in the constant, the treatment of each nonlinearity is slightly different and we simply present the details 
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for the readers' benefit. The similar error estimates of these four a-models show that, essentially, these 
four regularization models are consistent in their convergence to the exact 2D NSE system. 

Before we proceed in showing the main results of convergence rates of the a-modcls, we state an impor- 
tant lemma and a proposition which will play a crucial role in estimating the convergence rates of the 
various approximation models. 



\{ip-{I + a^Ay'cp,5)\ < 



Lemma 10. Let a > be a given fixed parameter, and let ip £ H and S £ V . Then 

l( 

Proof. First we observe that 

(p-{I + a^A)-^ip = a^A{I + a^A)-^(p. 

Therefore, 

|(^ - (/ + a^^)- V, '5)1 = \{a'A{I + a'A)-'^, 6)\ 

= \a (^{a^A)'/\l + a^A)-^^,A'/^6) \ 

< a\{a^Ay/'\l + a^Ay^ip\\\S\\ 

< a\\{a'A)^/^I + a'A)-^\\ciHM¥\\ 

<^M\si 

where the last inequality follows from the fact that the operator norm 

\\{a'Ay/^I + a'Ar^cm = (^^^) 

fc=l,2,... V 1 + " -^fc / 

y 1 

o<y<ooi + y 2 

which concludes our proof. □ 

Proposition 11. Let uq e D{A) and T > 0. Assume that a is small enough such that > 1, and 
suppose that is the solution of system (26) or (48) , then 

^sup^||u"(i)||ioe < ckl {l + \og (^)) . (51) 

Proof. By the Brezis-Gallouet inequality in Proposition 3, we have 

L \Au"'{t)[ 



-mu < c\K{t)r{i+\og 



27r ||u«(t)|| 



< c^iKwr + iKwriog^— j + iKwii^iog^ 

Therefore, by Proposition 6 for system (26) or Proposition 9 for system (48), we have a|AM"(t)| < Kq 
and ||u"(t)|| < I<o, for aU i G [0,T]. Since 3^ > 1, we obtain 

iK(^)iii^<ck^+i^o^iogr^U/?o^r^Vog^ 



Now, since by Proposition 6 or Proposition 9, ||u"(t)|| < A'o, for all t e [0, T], then we have 



A-„ J ° \\\u'{t)\\ - c- 
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Consequently, from all the above we have 



^sup^||«"(t)||i. < ckl {l + log 



□ 



Remark 6. The result above also holds for the solutions u" of systems (29) and (39) with A'o replaced 
by Kff2 in the cases of the NS-a and the ML-a models. 

Now we are ready to present the convergence rates of various a turbulence models. 

4.1. The rate of convergence of the Leray-a model. In this subsection, we will establish error 
estimates concerning the rate of convergence of the solutions of the 2D Leray-a model to solutions of the 
2D NSE, as the regularization parameter a goes to zero. We will proceed by estimating the L^-norm of 
the difference (5" = u — u", where u is the solution of 2D NSE system (22) and u" is the solution of 2D 
Leray-a model (26). 

From (22) and (26) we observe that S" ~ u ~ u" satisfies the following equation: 

—5" + lyAS" + B{u, (5") + B{5°',u) - B((5", (5") - (/ + A)-^ B{u°' .v"^) + B(w", u") ^ f - (I + a^A)-^ f. 
at 

(52) 

Theorem 12. Let u" be a solution of the 2D Leray-a system (26) with initial data ug S D{A), and let u 
he the solution of the 2D NSE system (22) with the same initial data uq over the interval [0,T]. Assume 
that a is small enough such that > 1 and let ^ u — , then 

sup |5"(t)P < e^ 

0<t<T 

where := (TKq (l + log (2:;;^)) + ll/lli,2((o t)-H) ) ' '^^^ ^'^^ given in Proposition 6 and 

5 respectively, and here c is a dimensionless constant that is independent of v, a, f or T . 

Proof. We take the inner product of the above equation (52) with 5" and use the identity (16) to obtain 

\J^\n'' + + iB{S^,u), 6") + ((/ + a'A)-\Biu'^,u'^) - B(u^,v^)), 6^) 

+ ((/ - (/ + a'A)'')B{u'^,u^),S^) = ((/ - (/ + a'Ar')f, 5^). 
Consequently, one has 

+ An? <\Ji\ + \J2\ + \h\ + (53) 

where 

Ji = (S(5",^.),<S"), (54) 
J2= ((/ + a'^)"'(S(?/°,w")-B(u°,i;")),,5"), (55) 
^3= ((/-(/ + a' A)-i)B(u", «"),,?"), (56) 
J4= ((/-(/ + a2A)-i)/,<5"). (57) 

Next, we estimate | Ji|, \ J2\1 l^al and | J4I. Thanks to estimate (11) and Young's inequality we have 

|Ji| = \{B{5^,u),5'^)\ < cnsm\u\\ < J||d-f + ^|5"p||^f . (58) 



16 



Y.CAO AND E.S. TITI 



Since it" + a^Au'^ = v", we obtain from (55) 

IJ2I = |((/ + a2A)-iB(u",a2Au"),<5")| 
= |(5(u",a2Aw"),(/ + a^A)^M")|, 

and by applying (15), (12) and Young's inequality we have 

IJ2I = |(B(w",(/ + a2^)-M"),a2Aw")| 

< ca2||w"||ioo||(/ + a2A)-i(5"|||Au"| 

< ca2||u"||L-|Au"|P"|| 

From Proposition 6, we have sup a^|^u"|^ < Kq, then by further applying Proposition 11 we have 

0<t<T 



IJ2I < 



a ||2 



9 

ca ~ , 



A'o* 1 + log 



\2TTa J 



(59) 



Now, we are ready to estimate | J3 1 . 

By Lemma 10, (56) and Young's inequality we have 



IJ3I < ^\Biu",u'')\\\6''\\<'^\\S"\\' + ^\Biu",u")\^ 

2v 



< JiFf + — !Klli==IKII^ 

o ly 



Since sup ||u"(i)|p < Kq by Proposition 6, then by Proposition 11 we reach 

Q<t<T 

(—) 

\2na J 



k3i< oP"r + — i^oMi+iog 



Finally, we estimate \Ji\. By Lemma 10, (57) and Young's inequality we obtain 



(60) 



\.h\<^\f\\\s'^\\<^\\n\' + ^\f\'- 

Now from (53), (58), (59), (60) and (61) we have 

V .iini^ < i|i.f r p + (i/p + A- (1 + log (^ 



(81) 



Dropping the i/||(5"|| term from the left-hand side and applying Gronwall's inequality we obtain 



|r(t)|2 < e^'/'o ll"('")ll''''"|(5"(0)|^ 



ca 

V 



2 ft 



|/(a)|2+i^„Ml + log 



— ) 

2-KOL J 



da. 



Since (5"(0) = 0, then by Propositon 5 we obtain 



\S'mf <—e^ [TK'Jl + log 



\L^{{0,T);H} 



which concludes our proof. 



□ 
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4.2. The rate of convergence of the Navier-Stokes-a model. In this subsection, we will establish 
error estimates concerning the rate of convergence of solutions of the 2D NS-a regularization model to 
solutions of the 2D NSE. as the regularization parameter a goes to zero. Wc will proceed by estimating 
the L^-norm of the difference 6" = u — u", where u is the solution of 2D NSE system (22) and is the 
solution of 2D NS-a model (29). 

From (22) and (29) we observe that 6" = u — u" satisfies the following equation: 

^5" + i^AS" + BiS", u) + B{u, S") ~ B{S°',5°') + {I + a'^A)-\B{u°',u°') ~ B{u",v")) 

+ (/ - (/ + a2yl)-i)s(u", u") + (/ + a^A)-^B*{u'^,v'^) = {I - {I + a^A)-^)f, (62) 

where for every v £V fixed, the operator B*{-,v) is the adjoint operator of the operator B{-,v), which 
is defined by {B* {ip,v), (p) = {B{f,v),ilj) for every <p,-i/> G V (see relation (18)). 

Theorem 13. Let u" be the solution of the 2D NS-a system (29) with initial data uq G D(A), and let u 
be the solution of the 2D NSE system (22) with the same initial data uo over the interval [0,T]. Assume 
that a is small enough such that > 1 and let 5°' ^ u — , then 

sup \S"it)f < 

0<t<T 

where ?2 ^e^ (tK*^ (l + log {^)) + ||/lli2((o,T);H)) ' ^nd Kq, Ko2 are given m (24) and (32), 
respectively. 

Proof. By taking inner product of (62) with 5" and using (16) and (17) we have 

i^|<5"P + + Jl + J2 + .h + J5 = J4, 

where Ji, J2,^3 and J4 arc given in (54)-(57) and J5 = ((/ + a^A)~^i?* (u", t;"), i5"). Therefore, 

+ An? <\Ji\ + 1 J2I + 1 J3I + + \M. (63) 

The estimates for | Ji|, | J2I, l^al and | J4I follow exactly the same steps as in the proof of Theorem 12, 
where Kq is replaced by 7^02 from Proposition 7 (see also Remark 6). 

Next, we estimate jJaj. To this end wc first observe that B*{u,u) = 0, because, thanks to (16), 
{B*{u, u),w)v' — {B{w, u),u)v' — for every w & V. 

Thus, 

IJ5I = |((/ + a2A)-ii?*(ri",«°),5")| 

= |((/ + a^A)-\B*{u", «") - B*(u", u")), <5")| 
= a^\{{I + a^A)-^B*{u",Au°'),S°')\, 

where we used the relation + a^Au" = v" and the bilinearity of B*{-, •). 
Consequently, by the definition of B* and (15), wc have 

I J5I = a^\iB{{I + a^A)-^S"',Au°'), «")[ a^\iB{{I + a^A)-!^", u"), Au°')\. 
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Therefore, by Lemma 2 and Young's inequality we have 

IJ5I < ca2||(/ + a2A)-i,5"||||u"|||Aw"| 
< ca^\\6"\\\\u''\\\Au"\ 

8 V 

By the above and Proposition 7 we conclude 

\J5\ < + —K^2- (64) 

From (58)-(61) (where Kq is replaced by Ko2; see also Remark 6), (63) and (64) we obtain 

^^\S-\^ + .WS'^r < + (l/P + (1 + log (^) 

Dropping the term from the left-hand side and applying Gronwall's inequality we obtain 

\S^{t)\^ < fo rf-|<5o(o)|2 + ^J^e^ ""^^^H'''^ (\f{<j)f + K^, (^1 + log (^))) da. 

Since (5"(0) = 0, then by Proposition 5 and the above we obtain 

,,2 .R-: 



r{t)\-' < (^tR'o2 (1 + log (^)) + WfWhmTyM) ) 



□ 



4.3. The rate of convergence of the Modified Leray-a model. In this subsection, we will establish 
error estimates concerning the rate of convergence of solutions of the 2D ML-a regularization model to 
solutions of the 2D NSE, as the regularization parameter a goes to zero. Wc will proceed by estimating 
the L^-norm of the difference S" = u — u", where u is the solution of 2D NSE system (22) and u" is the 
solution of the 2D ML-a model (39). 

From (22) and (39) we observe that (5" = m — u" satisfies the following equation: 

4-5" + B(r,u) + B(ii,(5")-B(5",r) + (/ + a'^)-i(B(w",u")-B(w",M")) 
dt 

+ (I-{I + a^A)-')B{u'',u'') = {I-{I + a'A)-')f. (65) 

Theorem 14. Let u" be the solution of the 2D Modified Leray-a system (39) with initial data uq G D{A), 
and let u he the solution of the 2D NSE system (22) with the same initial data uq over the interval [0,T]. 
Assume that a is small enough such that > 1 and let — u — u" , then 

sup |r(t)p < e^ 

0<t<T 

where e is given in Theorem 13. 

Proof. Taking the inner product of the above equation (65) with 5°' we obtain 

+ An? <\Ji\ + \h\ + I./4I + I Jel, (66) 
where Ji, J3 and J4 are as in (54), (56) and (57) respectively and 

Je = ((/ + a2A)-i(S(u",u") -S(z;",u")),,5"). 
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Here again the estimates for Ji, J3 and J4 are as in the proof of Theorem 12, with Kq replaced by K02 
(see also Remark 6). Next, we estimate Jq. Thanks to (14) we have 

\Je\ < |a2((/ + a^^)~^S(AM",M"),r)| 

= |a2(B(AM",(/ + a2A)-M"),M")| 

< ca^\\u"\\L^ + a^A)-^5°'\\\Au°'\ 

< ca2||it"||ioc|^u"|||,5"|| 

8 

Now since a^|ylu"p < Kq2 by Proposition 8, then by Proposition 11 (sec also Remark 6) we have 



I' , 



\-M<^^\\Sr + -K'o.[l + lo,^^)). (67) 
From (58), (60), (61) (where Kq is replaced by Kq2), (67) and (66) we have 

^m' + Hi*"f < ^ii""f i^-p + [up + 4'. (1 + 1=« (5I;))) . 

which concludes the proof of the theorem by Gronwall's inequality. □ 

4.4. The rate of convergence of the simplified Bardina model. Now we present the error esti- 
mates for the simplified Bardina rcgularization model. As before, we will proceed by estimating the 
L^-norm of the difference ^ u — u", where u is the solution of 2D NSE system (22) and u" is the 
solution of the simplified Bardina system (48). 

From (22) and (48) we observe that 5" = m — u" satisfies the following equation: 

^S°' + i^AS°' + B(u,6'') + B(6",u)^B(S°',S'') + [I - (I + A)-^)B(u°' ,u°') 
dt 

= {I-{I + a'A)-^)f. (68) 

Theorem 15. Let be the solution of the 2D simplified Bardina system (4-8) with the initial date 
uq e D{A), and let u he the solution of the 2D NSE system (22) with the same initial data uq over the 
interval [0,T]. Assume that a is small enough such that > 1 and let ^ u — u" , then 

sup 1^(01' <e^ 

0<t<T 

where e is given in Theorem 12. 

Proof. Taking the inner product of the equation (68) with 5" we obtain 

+ Hn? <\Ji\ + \m + \jaI (69) 

where Ji, J3 and J4 are as in (54), (56) and (57) respectively. Here again the estimates for Ji, J3 and J4 
are as in the proof of Theorem 12, with exactly the same Kq. 
From (58), (60), (61) and (69) we have 

j^\n + -ll^"f < ^ll-in-^-p + (l/P + A'o* (1 + log (^) 
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Dropping the term from the left-hand side and applying Gronwall's inequality we obtain 

Since (5"(0) = 0, then by Proposition 5 we obtain 



ca 



< [tkI U + log ) + ||/lli.,(o,T).m ) ^ 



V 



,2 



□ 

5. Error Estimates of The Galerkin Approximation of the a-RECULARizATioN Models 

In numerical simulation one needs to approximate the exact solutions of the underlying equations, 
that lie in the infinite dimensional function spaces, by functions that lie in finite dimensional spaces. In 
this section we estimate the errors between the approximation solutions u^, of the finite dimensional 
ordinary differential equation system (the finite dimensional Galerkin system of order m in this context) 
and the exact solutions of the a-models in L^-norm. The errors are given in terms of m and the 
regularization parameter a. Combining this with the results we establish in the previous section and by 
applying the triangle inequality, we obtain error estimates of numerical approximation solutions of 
the underlying a-model and of the exact solution u of the 2D NSE system. We take the Leray-a model 
as an example in this section and show the error estimate results. By similar arguments, one can show 
the error estimates for the other a-regularization models introduced in this paper. 

Galerkin approximation for the 2D Leray-a model 

Now we present the error estimates for the Galerkin approximation of the Leray-a model. This estimate 
for the rate of convergence of the Leray-a model is along the same lines of [17, 23, 38, 39] and in the 
spirit of the work for the 2D NSE. 

For an exact solution of the Leray-a system (26), we can decompose u" as follows: = p„j + 
where p„i = Pmu" and q„i = {I — Pm)u°' , Pm is the orthogonal project from H onto Hm, which is defined 
in Section 2. Rewriting equation (26) as 

— + vAu'^ + + a^A)-^B{u'^, {I + a^A)u") = (/ + a^A)-\f, (70) 

and since u" ~ p„i + gm, we can decompose the above equation (70) into the following coupled system 
of equations: 

^ + vAp„, + (/ + a2A)-ip,„ B(u«, (/ + a''A)u'^) = (I + a^Ar^P.^J, (71) 

% + vAq„, + (/ + a'A)-^! - P.„,)B{u^, {I + a^A)^") = (/ + a^A)-'{I - P^)/. (72) 
dt 

For the Galerkin approximation system of the Leray-a system, we rewrite the equation (27) in the 
following equivalent form 

+ vAuZ, + + a^A)-^P^B{uZ,, {I + a^A)u'^) = (/ + a^A)-^P„J. (73) 

We will proceed by first estimating the L^-norm of qm and then the L^-norm of the difference 5m = 
Pm — where is the solution of the Galerkin system of the Leray-a system (73). Then by the trian- 
gle inequality and orthogonality of spaces projected by Pm and (/ — -Pm), we obtain the error estimates 
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of the L^-norm: |u" - <J-^ = + \pm ~ <r 

From (71) and (73) we observe that 5m = Pm — satisfies the following equation: 

^ + vASra + {! + a^A)-^PmB{u^, {I + a^A)!^ - {I + a^A)-' P,nB{iC {I + a^A)u^) = 0. 
Since u" = p„i + Qm and Sm = Pm — u"^, we can rewrite the above equation as 
^ + ,^A6^ + iI + a^Ay^PmB{6m + qm,iI + c^^A)u") 

- (/ + a^A)-'PraB{u^, {I + a^A){Sm + g™)) - 0. (74) 

Theorem 16. Let T > and let u" be a solution of the Leray-a system (70) with initial data uq G D^A), 
and let u'^ be the solution of (73) with initial data UQm ~ Pm uq over the interval [0,T]. For a given 
TO > 1, assume that a is small such that o? < — , then 

— — A„, + i 

sup |M"(t)-<,(t)p <e2, 

0<t<T 



where 



■■= j^^iQ + R + LmUV), im = l + log(^) andQ,R,U,V depend on y,uo,\og{^), f ,T , 
and are given explicitly in (80), (81), (92) and (93), respectively. 

Remark 7. Here we require / e L°°{{0,T); H), which is stronger assumption than the condition / S 
L^((0, T); H) in the estimate of rates of convergence in section 4. 

Proof. First, we estimate the norm of qm- 

We take the inner product of equation (72) with qm and get 

— kmf + ^hmW < m + a^A)-'Biu'^,iI + a^A)u'^),qm)\ 

+ \{{I + c^^A)-\f,qm)\ 
< M1+M2. (75) 

Next, we estimat Mi, Af2- By virtue of (12) we have 

Ml = + a^A)-^B{u",{I + a^A)u''),qm)\ 

= |(B(u", (/ + a2A)u"), (/ + a^A)-\m)\ 

< c\\u"\\l^ + a^A)u''\\\{I + a^Ay'qml 

< cKoil + log(^))V2(||^"|| + «2|^3/2^"|)^ , 

where in the last inequality we apply (5) and Proposition 11 to ||w"||ioc. Notice that by the assumption 
that < x^'tT — XT ~ (^)^' consequently have > 1, and as a result, it is valid to apply 
Proposition 11. 

Now, by Young's inequality and the a priori estimates obtained in (28), we have 

Ml < J|<z™f + -^-A'o^(l + log(-^))(h"f + 
4 Am+i 27ra 

< ^|k™f + -^-/^o'(l + log(-^)) (76) 
4 A,„+i at: a 

-^i^„2(l + log(-^))(a2|^3/2."n. (77) 



Am+i ly 27ra 
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Now, for estimating Af2, we apply Cauchy-Schwarz inequality, (5) and Young's inequality to obtain 
M2 = |((/ + a2A)-i(/-P„)/,g„)| = |((/-P„)/,(/ + a2^)-ig„)| 

< |/||(/ + «2yl)-i,,„|<|/|" <^||,„||2+ 1 (78) 

Let us substitute the bounds for Afi, M2 into (75) to obtain 

' |2 , .,11 112 / 1 C / 2 , r>4^- ■ ' 



£At(l + log(-^))(a2|A3/2."n. (79) 



Am+l 27rQ! 

By Poincare inequality and (5), A„i+i|(7„ip < Ij^mP, we obtain 

+ T^-/^o'(l+log(^))(a^l^^/^-"n. 

By Gronwall inequality, we get 

|g™(t)P < e-^^'"+^*|g„(0)p+ Te-'^^^+^f'-^^Y^- f|/P + At(l + log(-^-))) ds 

Jo Am+i V \ 27ra / 

*e-''^'"+i(*-^'-^-At(l + log(— ))(a2|A3/2u"|2) 

Arn+l V 27ra 



< e-W.*LlgiM + T2^(l - e-''^™-^^ (l/li=^((o.T);^) + ^0^(1 + log(^)) 

^m+l ^m+1 V 

+ -^^A:o'(l+log(-^)) Aa^i^^/V^d. 
Am+l V lira Jq 

< -^lAuol' + T5^-J (\f\Uio,Ty.H) + + log(^)) 
'ri^o^(l + log(;^)). 



Am+l ' '^^'•2TTa' 

where we apply the a priori estimates of solutions of Leray-a model given in Proposition 6. Next, denote 

by 



Q \Auo\' + ^ (|/li»((o,T)://) + ^0(1 + log(^)) ) ' (80) 



and 



we obtain 



i?:=__At(l + log(^)). (81) 



-^m+1 Am+l 



By the assumption that < x^^j we have 

1 

'^m+l 



9ml' <T2 (Q + R). (82) 

A, 
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Next, we estimate the L^-norm of Sm- 

Taking the inner product of equation (74) with Sm, we get 



^^|<5™P + i^\\5m\\^ < Ms + A'h + M5 + Me, (83) 



where 



M3 = \Bi6rn, (/ + a^A)u''), {I + a^A)-'drn)\, 
M4 = \B{q,n, (/ + a^A)?/"), (/ + a^A)~Hm)\, 
Ms = \{B{u^, {I + a^A)6„,), {I + a^Ay^6„,)\, 
Me = \{B{u^, {I + a^A)qm), {I + a^Ay^S^)]- 

First, let us estimate Thanks to (11) we have 

A/3 < c|<5™|i/2p„||i/2||(j + a'A)u''\\\iI + a2A)-i5„|i/2||(/ + a'Ay'SmW'^' 

< c|<5™|||5™||(||«"||+a2|A3/2w-|) 

< ^\\Smr + -\S,n\'i\\u-r + a^\A'/'u-\') 
16 

< ^\\S,nr + -\S,n\'K'o + ■^-ISmWa'lA'/'u-l^), (84) 

16 \,n+l 

where in the last inequahty we use the a priori estimates of Leray-a obtained in Proposition 6 and the 
assumption that < ^ - 

Now, we estimate M4. By applying (21) and (7) we obtain 



Af4 < c\qm\\\iI + a^A)u''\\\\iI + a^A)-'Sm\\i 



< c\qm\i\\u^\\+a^\A^/'u"\)\\{I + a^A)-'5m\\ [l + \og{^) 

< c\qm\{\K\\+a'\A'/'u-\)\\S^\\Ll/^ 

< ^\\Sr.r + —iqMww^r + «^(a^iA3/2^"n) 

16 ly 
16 ly 

where Lm = 1 + log(^), and we apply Proposition 3 and (7) in the second inequality above. 
Let us now estimate By virtue of (11) we have 

A/5 < c|7.j;|i/2||^- ||i/2||(/ + a'A)Srn\\\{I + A) " | || ^ a2A)-M™||i/2 
< cK|i/2|K||i/2(l + a2A™)|||<5„|||5™|i/2||5„||i/2. 

Since a^Xm < ct^Am+i < 1, and c denotes general dimensionless constant, we have 

A/5 < c||<5,„|l3/2|<5,„|i/2|^^ji/2|j^.|ji/2 

16 ly 

< + (86) 
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where in the last inequahty we apply Proposition 6. 

For the last term Afg, we will proceed by applying (15) and (12), then we obtain 

Me = |(B(M^,(/ + aM)g,„),(/ + a2A)-i,5™)| 
= |(B«„ (/ + a^A)-^6„,), {I + a^A)q,„,)\ 
< cIKJIloo ||(/ + a^A)~'S,n\\\{I + a^A)q.^\. 
By applying Proposition 3, inequality (7), (5), and the facts that 

\\{I + a'^A)-^4)\\<\\(f)\\, ioraWcfxEV, 

\{I + a^A)4>\ - \(j)\+a^\A(j)\ for aU e D{A), 



we obtain 



Me < clKJI (l + log(^)j \\{I + a'A)-^S,J\iI + a^A)qm\ 



< cLll'^IU 1 1 (5™ 1 1 ( I I + | Aq„ 

< cLl^^K,\\S,n\\{\qm\+a^^-^^ 



1/2 
m+1 



16 A„i+i 

< ^||<5,„|r + ^i-|g„|^ + -^^/^o^(aV^/^-"^), (87) 

where i„i = 1 + log(^), and in the second inequality above we use the fact that |[it"J[ < A'o from 
Proposition 6. 

Plugging (84), (85), (86) and (87) into the inequality (83) we obtain 

^M' + .\\5,nf < ^|5™|^(A'o' + ^(aV^/^-"r) + 4f) 

at V -^m+l V A\ 

cLm £>2| |2 I |2/ 2|^3/2 Q!|2\ 

+ K^lq^l +a \qm\ [a lA-"' u ) 

1 cLr. 



|2 



/?2(a2|^3/2^a|2) 



where 



and 



V{t) -.^ £a^2|^^|2^^2C|^^^|2(^2|^3/2^.|2)^ 1 £^2(^2|^3/2^a|2)^ (gg) 



Applying Gronwall inequality and recalling that |(5,„(0)| = 0, we obtain 



By the a priori estimates in Proposition 6, we know that 

V / (a2|A3/2w"|2)dt < (91) 
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which imphes that U G L^{0,T) and consequently we have 

r U{t) dt < -{klT + -^kl + =■■ U. (92) 

Now for estimating Vit)^ we recall from (82) that \qm\^ < -vr- — {Q + R)- Applying (91), we obtain 

f V{t) dt < -J— {-kliQ + R)T + -^^kliQ + R) + -^ki] =: -^V . (93) 
Plugging back into the inequality (90), we get 



\5nm < TT^UV. (94) 

\n+l 



Now by Pythagorean Theorem and orthogonality of the projection spaces P,„ and (/ — Pm), we have 

1 

0<t<T Q<t<T 0<t<T 

where Q,R and U,V are given in (80), (81), (92) and (93), respectively. □ 



sup |w"(t) -<(t)|" < sup \Srn\''+ sup \qrnf<^^iQ + R+Lm.UV), (95) 



Remark 8. The above result |w"(t)— w"j(t)| = O (^ ^ "^l^ j is of the same order as that of the 

error estimates for the usual Galerkin approximation of NSE. Indeed, in [44] the author points out that 
for the 2D NSE, the error estimate for the usual Galerkin, \u~Um\ is of order O ^ j^^^(log(Ami^))^/^)^ , 

where u is the solution of 2D NSE (22) and u,n is the solution of the corresponding Galerkin system (23). 
Furthermore, this estimate is optimal up to the logarithmic terms provided f G L^. 



Now, simply by applying the triangle inequality, we achieve the error estimates of the solution for 
the finite-dimensional Galerkin system of the 2D Leray-a model as approximation of the exact solution 
u of the 2D NSE system. 

Theorem 17. Let T > and let u he the solution of the 2D NSE system (22) with initial data uq € D{A) 
and be the solution of finite- dimensional Galerkin approximation of the 2D Leray-a system (73) with 
initial data u^m = Pm ua over the interval [0, T] . For a given m > 1, assume that a is small enough such 
that a < Y"^f7L ' ^'^s^s Kn+i is the (m + l)-th eigenvalue of the Stokes operator A, then 

sup Kf)-^^(Op<cf(-^^)Mogf^^^)), (96) 
o<t<T \ \m+iL^ \ (27r)^ // 

where C is a constant, which depends on uq, f, L and T only. 

Proof. By the triangle inequality we have 

\u{t) - <,(0P < 2(1^.(0 ~ u-{t)\'' + \u-{t) - <,(t)p). (97) 

From Theorem 12 we obtain, under the assumption that a < 

.S,i»("-"°('>i'^'''((^)'"'8(2^;j))' w 

where C" is a constant that depends only on uq, f, L and T. By virtue of Theorem 16 above, we have, 
under the assumption that < — , 



+1 




^sup^|h"(0 - u^(t)|^ <C{[ -^-^ ) log ( ) ) , (99) 
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where C is a constant that depends only on uo, /, L and T. 

Now we assume that a < x^frx- Noticing that Ai = {^)^, we have 

, , 27r 27r , Ai 1 1 , , 

" - ^A T^^A ^ A - A ' 

where in the last step we use the fact that Ai < Xm+i- 

Since both assumptions for (98) and (99) are satisfied, we combine (97), (98) and (99) to obtain 
sup |u(t) -u^(t)p < C 

0<t<T 

which concludes our proof. □ 
Recall that, according to (1), the (m + l)-th eigenvalue of the 2D Stokes operator satisfies 

m + 1 A,„+i , , , -,^ 

<— — <co(m + f), (101) 

Co Ai 

where cq is a constant that depends only on L. Applying the above asymptotic estimate and the fact 
Ai = (x")^' '-"^^ '^^'^ rewrite the conditions and results of Theorem 17 in the form of the following 
corollary: 

Corollary 18. Let T > and let u he the solution of the 2D NSE system (22) with initial data uq G D{A) 
and be the solution of finite- dimensional Galerkin approximation of the 2D Leray-a system (73) with 
initial data uom = Pm uo over the interval [0, T] . For a given m > I, assume that a is small enough such 
that a < ^ ^^^^-^ , where Cq is a constant that depends only on L as stated in (1), then 

sup \u[t)~uZAt)? < Co(^-)2log(m+ 1) (102) 
Q<t<T m + 1 

where Cq is a constant, which depends on u, uq, f, L and T only. 

The error estimates for the other a-models arc not presented in this paper in detail. The results are 
the same as that of the Leray-a model we present above. One can follow the same idea and check the 
calculation. In fact, the proof of the estimates for the Modified Leray-a model and simplified Bardina 
model will follow readily from the proof of the Leray-a model since the former two have milder nonlinearity 
and the estimates on the nonlinear terms will be easier compared to that of the Lcray-a model. For the 
2D Navier-Stokes-a equation, combining (12)-(14) and (18), we can get similar inequalities for B as 
those of B and by following same steps obtain similar error estimates for the Navier-Stokes-a model. In 
a subsequent paper we will present error estimates regarding the three-dimensional case in the spirit of 
the results reported in [9] . 

APPENDIX 

For the sake of completeness, we present in this section the version of the Brezis-Gallouet inequality 
presented in Proposition 3. 

Proof. Let if G D{A) such that 

/cGZ2\(0,0) 



therefore, 



fcGZ=^\(0,0) 0<|fc|<A/ |fe|>A/ 

^ + , , ^ (i + ifci 



2^ 



0<1M<A/ (1 + |/fc|2)l/2 I'^'^l + |fc|t:M (1 + |fc|2) 1^'=! 
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where M > 1 is a real number to be determined later. 
Then by the Cauchy-Schwarz inequality we have 

1/2 / 1 \ 1/2 



0<|fe|<A/ / \0<|fc|<A/ (1 + |fcp) 



1/2 / ^ \ 1/2 



|fc|>J\/'' ' ' J \\k\>M {1 + Ikl"^)^ 



\ 1/2 / V 1/2 



I\V\<M 

We observe that in two dimensions 



and 



Therefore, 



J\y\<M^ + \yr Jo Jo 1 + 7-^ 



dy Z"^" ,^ f°° rdr 1 



II^IU. < c (m (log(l + A/^))^/% 



for every M > 1. Notice that by Poincare inequality we have > 1, therefore, one can choose 

M = + 1 to conclude the proof. □ 
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